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SUMMARY 


\ 

\ 

\ 

\ 

We  derived  closed-form  expressions  for  the  Cramer-Rao  Bound,  MUSIC,  and  Maximum 
Likelihood  (ML)  asymptotic  variances  corresponding  to  the  two-source  direction-of-arrivaj 
estimation  where  sources  were  modeled  as  deterministic  signals  impinging  on  a  uniform 
linear  array.  The  choice  of  the  center  of  the  array  as  the  coordinate  reference  resulted  in 
compact  expressions  that  greatly  facilitated  our  study  of  effects  of  temporal  phase  differ¬ 
ence  (correlation  phase)  of  the  two  sources  on  asymptotic  variances  of  estimation  error. 
These  effects  were  shown  to  be  intensified  when  the  two  signals  were  closely  spaced  and/or 
when  their  normalized  correlation  magnitude  was  high.  Numerical  examples  obtained  from 
specializing  general  formulas  agreed  with  our  results.  /  yn 
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KEY  SYMBOLS 


m  =  number  of  sensors 

n  =  number  of  arrivals 

N  =  number  of  snapshots 

aj(u>j),  1  <  j  <n  :  m  x  1  steering  vectors 

dj  =  da.j(ujj)/dojj,  1  <  jf  <  n 

n  =  (u>i,uj2,  ...  ,un)T  :  parameter  vector 

A  =  ai(u>i):  •  •  •  :’an(cjn  )|  :  m  x  n  transfer  matrix  (direction  matrix) 

D  =  dii  •  •  •  :dn 

A+  =  (A*A)-1A*  :  left  pseudo-inverse  of  A 

PA  =  AA+  :  orthogonal  projection  on  subspace  spanned  by  columns  of  A 

Pa  =  1  ~  pa 
H  =  D*P^D 

S  =  n  x  n  source  covariance  matrix 


R  =  ASA*  +  o2I  =  ESA,E*  +  E„A„E*  :  m  x  m  covariance  matrix  of  array  output 


>  •••>  Xd>  -Wi  =  •  •  •  =  Am  =  o'2,  d  <  n  :  eigenvalues  of  R 
E*  =  jei :  •  •  •  iejj  :  columns  of  E,  are  normalized  eigenvectors  corresponding  to  the  d 
largest  eigenvalues  of  R 

E n  =  jej+1:  •  •  •  :em  :  columns  of  En  are  normalized  eigenvectors  corresponding  to 
the  eigenvalue  <x2  of  multiplicity  m  -  d 


M~r  =  (M-1)r  =  (Mt)-1 


M  0  N  :  Schur  (Hadamard,  elementwise)  product,  (M  0  N)tJ  =  M,_,Nt'j 
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1.  INTRODUCTION 


It  is  well  recognized  (references  [1  j ,  [2],  [3],  [4],  [5],  [6],  [7])  that,  for  many  important 
signal  processing  applications,  the  pertinent  problem  is  estimation  of  the  parameters  of  the 
following  basic  model: 

x(t)  =  As(<)  +  n(t)  (1) 

where  x(t)  is  an  observed  data  vector  of  size  m,  with  m  equal  to  the  number  of  sensors, 
the  n-vector  s (t)  contains  the  complex  envelopes  of  the  n  narrowband  signals  from  far-held 
emitters,  n(<)  denotes  a  complex  m-vector  of  additive  noise,  and  the  columns  of  the  m  X  n 
transfer  (direction)  matrix  A(Q)  are  steering  vectors  aj(uj),  1  <  j  <  n,  with  the  unknown 
parameter  vector  fl  =  (uq,^, . . .  ,u>n)T  corresponding  to  the  unknown  directions  of  arrivals 
(DOAs). 

The  waveforms  s(t)  are  modeled  as  stationary,  jointly  Gaussian  processes  with  zero 
mean  and  second  moments 

E  (s(t)s*(s))  =  S 6t,„  E  (s(t)sT(s))  =  0, 

where,  as  usual,  the  operator  E(-)  denotes  the  expected  value,  the  superscript  T  for  trans¬ 
pose,  the  superscript  *  for  complex  conjugate  and  Stf3  is  the  Kronecker  delta. 

Similarly,  the  additive  noise  n(<),  uncorrelated  with  the  signal  waveforms,  is  assumed 
to  be  a  stationary,  zero  mean,  Gaussian  process  with  second  moments 

E(n(t)n*(s))  =  <r2I6t,3 ,  E  (n(t)nr(s))  =  0. 

The  array  output,  being  observed  at  N  discrete  time  instances,  is  a  stationary,  zero 
mean  Gaussian  process  with  second  moments 

E  (x(t)x*(s))  =  R St,,  =  { A(Il)S  A*(n)  +  E  (x(/)xr(s))  =  0. 

It  is  well  known  that  R  has  the  following  eigen-decomposition 

m 

R  =  Aje/e*  ~  E,AjE‘  +  EnAnE*, 

i=i 

with 

Ai  A  •  •  *  ^  Xd  }>  A^j  —  *  —  Am  —  ^  4  7i, 

where  the  columns  of  E*  are  the  eigenvectors  of  the  corresponding  d  largest  eigenvalues  and 
the  columns  of  En  are  the  n-d  eigenvectors  of  the  eigenvalue  a 2.  Here,  the  diagonal  matrix 
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A,  contains  the  d  largest,  in  descending  order,  eigenvalues  of  R.  The  range  of  E3  is  often 
referred  to  as  the  signal  subspace  and  its  orthogonal  complement,  which  is  spanned  by  En, 
as  the  noise  subspace. 

The  problem  of  determination  of  the  number  of  arrivals  n  belongs  to  the  area  of  detection 
and  will  not  be  discussed  here.  It  is  assumed  that  n  and  d  are  known  and  one  wishes  to 
estimate  the  DOA  vector  SI. 

Recently  there  has  been  interest  in  treating  the  matrix  S  other  than  real  for  two-source 
applications  (references  [8],  [9],  [10],  [11],  [12]).  Earlier  analytical  results  in  this  direction 
for  the  Cramer-Rao  bound  can  be  found  in  references  [13],  [14],  [15],  and  most  recently  in 
reference  [16].  The  purpose  of  this  report  is  to  investigate  effects  of  temporal  phase  difference 
(that  is,  correlation  phase)  of  two  sources  on  asymptotic  variances  of  estimation  error  using 
a  uniform  linear  array.  In  particular,  we  shall  restrict  to  the  Cramer-Rao  bound  (CRB), 
the  MUSIC  and  Maximum  Likelihood  (ML)  methods  that  are  of  practical  and  theorectical 
interest  (references  [2],  [3],  [4]).  Our  results  also  provide  an  analytical  background  for 
simulations  such  as  those  performed  in  references  [10]  and  [12]. 

This  paper  is  organized  as  follows.  In  section  2,  we  shall  collect  several  important 
asymptotic  results  on  estimation  error  for  the  deterministic-signal  model,  i.e.,  in  (1)  the 
waveforms  s(t)  are  fixed  in  all  realizations  of  the  random  data  x(t).  Among  these  are  the 
CRB  for  the  estimation  of  the  parameter  vector  f 2  of  (1),  the  asymptotic  covariance  matrix 
for  the  MUSIC  method  and  that  for  the  ML  method  -  all  three  were  derived  by  Stoica 
and  Nehorai  (references  [2],  [3]).  Recently,  the  asymptotic  covariance  matrix  for  a  general 
multidimensional  signal  subspace  method,  which  includes  ML,  has  been  obtained  by  Viberg 
and  Ottersten  (references  [5],  [6],  [7]).  For  an  appropriate  choice  of  the  weighting  matrix, 
the  resulting  covariance  matrix  of  estimation  error  is  smallest  among  those  in  the  weighted 
subspace  fitting  (WSF). 

We  shall  specify  the  geometry  of  the  sensor  array  and  the  number  of  arrivals  in  section 
3.  Precisely,  we  shall  study  the  situation  of  two  narrowband  plane  waves  impinging  on  a 
uniform  linear  array  (ULA).  This  will  specify  the  form  of  the  steering  vectors,  and  thus 
the  transfer  matrix.  Basic  properties  of  ULA  that  are  essential  for  subsequent  development 
are  derived.  Some  of  these  properties  are  independent  of  coordinate  reference;  others  are 
not,  and  an  inappropriate  choice  will  result  in  unnecessarily  complicated  expressions  and 
analysis.  A  suitable  coordinate  reference  for  our  purpose  is  the  center  of  the  array. 
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In  section  4,  we  shall  be  concerned  with  effects  of  the  temporal  phase  difference  on  the 
CRB  covariance.  The  main  result  in  this  section  states  that,  for  the  coordinate  reference  at 
the  center  of  the  array,  then  the  CRB  variance  of  each  estimation,  considered  as  a  function 
of  the  phase  difference  0,  is  periodic  of  period  180°,  symmetric  about  90°,  and  on  the  interval 
where  0°  <  0  <  180°,  it  is  decreasing  when  0°  <  9  <  90°,  increasing  when  90°  <  6  <  180°, 
and  it  assumes  the  minimum  at  90°  and  maximum  (either  finite  or  infinite)  at  0°  and  180° 
(that  is,  either  the  two  signals  are  in  phase  or  out  of  phase).  This  behavior,  observed  in 
simulation  studies  (references  [10],  [12]),  can  be  malytically  deduced  from  reference  [15]. 
Our  proof  follows  the  development  in  reference  [2]. 

Results  on  effects  of  the  temporal  phase  difference  on  asymptotic  MUSIC  covariance 
is  derived  in  section  5.  We  analytically  show  that,  for  the  coordinate  reference  at  the 
center  of  the  array,  the  asymptotic  MUSIC  variance  of  each  estimation  is  of  the  form 
A  +  B  cos(0),  where  A  >  |i?|  >  0  with  B  having  the  sign  of  the  function  p(Aw;m)  = 
••  (The  electrical  angular  separation  Au  will  be  defined  subsequently  in  section  2 
below).  The  behavior  of  the  asymptotic  MUSIC  variance  with  respect  to  9  thus  can  be 
obtained  immediately.  In  particular,  suppose  the  electrical  angular  separation  is  within  one 
beamwidth  (0  <  Aw  <  —),  then  the  asymptotic  MUSIC  variance  is  symmetric  about  180°, 
decreasing  when  0°  <  8  <  180°,  increasing  when  180°  <  6  <  360°,  attains  its  maximum  at 
9  =  0°  (that  is,  the  two  signals  are  in  phase)  and  minimum  at  9  =  180°  (that  is,  the  two 
signals  are  out  of  phase). 

In  section  6,  we  examine  effects  of  the  temporal  phase  difference  on  the  asymptotic  ML 
variance,  which  is  derived  in  terms  of  CRB  variances  and  covariance.  As  in  the  case  of  the 
asymptotic  MUSIC  variance,  the  behavior  of  the  asymptotic  ML  variance,  as  a  function  of 
9,  is  also  influenced  by  Aw  in  terms  of  beamwidths.  In  particular,  suppose  the  electrical 
angular  separation  is  within  one  beamwidth  (0  <  Aw  <  |£),  then  the  asymptotic  ML 
variance  is  symmetric  about  180°,  and  on  the  interval  where  0°  <  9  <  180°,  it  attains  the 
maximum  at  9  =  0°  (that  is,  the  two  signals  are  in  phase)  and  minimum  at  9  =  90  with 
90°  <  9o  <  180°,  where  the  exact  location  of  8o  depends  on  Aw  and  other  parameters. 

Numerical  results  obtained  from  specializing  general  formulas  surveyed  in  section  2  are 
presented  in  section  7.  These  numerical  results  are  consistent  with  our  analytical  results. 
As  discussed  in  sections  4  to  6,  effects  of  temporal  phase  difference  on  CRB,  asymptotic 
MUSIC,  and  ML  variances  of  estimation  “r^nr  are  intensified  when  t>  two  sources  are 


3 


closely  spaced  and/or  when  their  normalized  correlation  magnitude  is  high.  This  behavior 
is  well  noticeable  in  figures  1(a)  through  3(d)  . 

Finally,  in  section  8,  we  conclude  the  paper  by  summarizing  our  findings  with  emphasis 
on  the  case  where  the  two  arrivals  are  closely  spaced  with  their  electrical  angular  separation 
within  one  beamwidth. 

2.  ASYMPTOTIC  COVARIANCE  MATRICES  OF 
ESTIMATION  ERROR 

The  Cramer-Rao  bound  (CRB)  provides  a  lower  bound  for  the  covariance  matrix  of  the 
estimation  error  of  any  unbiased  estimate.  Stoica  and  Nehorai  (reference  [4])  gave  a  general 
formula  for  the  CRB  from  which  they  derived  expressions  for  the  case  where  signals  are 
deterministic  as  well  as  the  case  where  signals  are  random.  Since  the  expression  of  CRB 
for  random  signals  is  much  more  complicated  than  its  deterministic  counterpart,  and  since 
deterministic  CRB  is  the  lower  bound  of  random  CRB  (reference  [4]),  hereafter  we  shall 
restrict  our  attention  to  the  deterministic-signal  model,  i.e.,  in  (1)  the  waveforms  s(t)  are 
fixed  in  all  realizations  of  the  random  data  x(t)  and  only  n(t)  varies  from  realization  to 
realization.  For  deterministic  signals,  Stoica  and  Nehorai  proved 

Theorem  1  (reference  [2])  Let  (is  be  an  (asymptotic)  unbiased  estimate 
parameter  vector  O.  For  large  N,  the  Cramer-Rao  inequality  can  be  written 

NE  (((is  -  n)((ls  ~  n)T)  >  CRB, 

where 

CRB  =  y  [Re(H0ST)]_1. 

Here  the  matrices  D  and  are  evaluated  at  S7. 

Hereafter,  we  shall  refer  to 

lim  NE  (((is  -  Sl)((ls  -  «)T)  (3) 

as  the  asymptotic  covariance  matrix  of  the  estimation  error  corresponding  to  the  method 
that  yields  the  estimate  (l  of  the  true  parameter  vector  ft.  The  matrix  Re(H  ©  ST)  is 
called  the  Fisher  information  matrix. 

The  asymptotic  covariance  matrix  for  the  MUSIC  method  was  also  obtained  by  Stoica 
and  Nehorai;  in  matrix  form  it  reads 


of  the  true 


(2) 
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Proposition  1  (references  [2],  [3])  Suppose  the  source  covariance  matrix  S  is  of  full  rank. 
Then  the  asymptotic  covariance  matrix  of  the  estimation  error  of  the  MUSIC  method  is 

CMU  =  yCHOir^HoK^HOl)-1,  (4) 

where 

K  =  S"1  +  ct2S-1(A*A)-,S-1.  (5) 


Notice  that  (H  ©  I)  1  is  a  fancy  expression  of  the  matrix 


l 

TUT 


l 

h22 


1 

HTT 

Recently,  Viberg  and  Ottersten  (references  [5],  [6],  and  [7])  have  introduced  a  unified  ap¬ 
proach  to  multidimensional  signal  subspace  method  called  weighted  subspace  fitting  ( vVSF) 
by  incorporating  a  general  Hermitian  weighting  matrix  W  and  the  estimation  error  covari¬ 
ance  is  minimized  with  respect  to  this  weighting.  When  W  =  I,  the  corresponding  method 
is  called  the  multidemensional  MUSIC  (MD-MUSIC),  as  opposed  to  the  conventional  (one- 
dimensional  search)  MUSIC. 

The  maximum  likelihood  (ML)  method  is  known  (references  [5],  [6])  to  be  equivalent  to 
the  weighted  subspace  fitting  method  with  the  weighting  matrix  W  =  A  where 

A|  -  <72 

A2  —  o2 


1  A 


A  =  A,  -  crl  = 

A d  -  <72 

The  asymptotic  covariance  matrix  of  the  estimation  error  of  the  subspace  fitting  method 
with  the  weighting  matrix  W  is  given  by 


Theorem  2  (reference  [6])  For  the  subspace  fitting  method  with  the  Hermitian  positive 
definite  weighting  matrix  W,  the  asymptotic  distribution  of  the  estimation  error  is  given  by 

-  ft)  ~  N(0,  Civ sf), 
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where 


(6) 


CW  =  (V")-1Q(V,/)-1 

with 

V"  =  -2Re{H0(A+E,WE;A+*)r},  (7) 

Q  =  2<r2Re{H  0  (A+E,WAsA-2WE*A+*)t}.  (8) 

Here,  all  expressions  are  evaluated  at  the  true  parameter  Vt.  The  result  holds  for  arbitrary 
signal  correlation  including  full  coherence. 

When  S  has  full  rank,  a  convenient  expression  for  C ml  can  be  derived  from  (6)  to  (8): 

Corollary  1  (reference  [6])  Assume  the  source  covariance  matrix  S  is  nonsingular.  Then 
the  asymptotic  covariance  matrix  of  the  estimation  error  of  the  ML  method  (W  =  A)  is 

CML  =  CRB  [i  +  2Re{H  ©  (A*A)~r}CRB|  .  (9) 

This  expression  is  consistent  with  an  earlier  expression  obtained  by  Stoica  and  Nehorai 
in  reference  [3]. 

3.  BASIC  PROPERTIES  OF  UNIFORM  LINEAR  ARRAY 

In  this  section,  we  shall  discuss  fundamental  properties  of  a  ULA  that  are  pertinent  to 
subsequent  development.  Some  properties  are  independent  of  the  choice  of  coordinate  refer¬ 
ence  (Lemmas  1  and  2);  others  are  not,  and  an  unsuitable  choice  will  result  in  unnecessarily 
complicated  expressions  and  analysis.  An  appropriate  coordinate  reference  for  our  purpose 
is  the  center  of  the  array,  which  will  become  clear  by  the  end  of  the  section. 

We  shall  consider  the  situation  of  two  narrowband  plane  signals  impinging  on  a  uniformly 
distributed  linear  array.  In  this  case,  the  parameter  vector  is  R  =  (wi,W2)t-  In  view  of 
reference  [17],  we  shall  assume  the  number  of  array  elements  m  >  3. 

Let  Au  =  jwi  -  u>2i  be  the  electrical  angular  separation  of  the  two  signals  that  will  be 
defined  shortly.  In  this  section,  we  shall  prove  that  the  asymptotic  covariance  matrix  CRB, 
Cmu,  C wsf  (and  thus  C ml  and  Civ0)  are  dependent  on  the  electrical  angular  separation 
Aw  and  independent  of  the  absolute  values  wj  and  u>2-  An  analogous  temporal  result  for 
CRB  was  shown  in  references  [13]  and  [14].  Proceeding  from  spatial  premises  of  the  basic 
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mod<’  stated  in  (1),  our  development  gives  a  direct  treatment  and  insight  to  the  two-source 
direction  finding  (DF)  using  a  uniform  linear  array.  Consider  the  general  case  where  the 
coordinate  origin  O  >s  a*  a  point  on  the  line  segment  connecting  the  first  sensor  and  the 
last  (m-th)  sensor.  In  this  setting,  if  the  coordinate  origin  Q  is  between  the  p-th  sensor 
and  the  ( p  -1-  l)-st  sensor,  for  some  1  <  p  <  m  -  1,  and  let 

a  distance  between  Q  and  p-th  element  ,  ,  .  . 

r  = - : - - — — - - ,  where  0  <  r  <  1,  (10) 

interelement  spacing 

then  the  steering  vector  aj(u>j),  for  j  =  1,2,  is  given  by 


ai(wi)  =  [< 


-i(p-l+r)u-j 


p-ill+rjo/j  -i  rioj 

C  ,  C  . 


Pi(l-r)^ 


where 


■  (A.  \ 

ojj  =  —  sin(d>j) 


denotes  the  electrical  phase  angle  corresponding  to  the  angle  of  arrival  (f)j,  measured  with 
respect  to  the  normal  to  the  array,  of  the  j-th  plane  wave  with  d  as  the  interelement  spacing 
and  A  as  the  array  wavelength. 

The  case  where  the  coordinate  center  is  chosen  to  be  at  a  particular  sensor  element  can 
be  obtained  from  the  general  formulation  by  setting  either  r  =  0  or  r  =  1.  A  convenient 
choice  of  coordinate  origin  of  the  array  is  its  “center,”  which  is  its  actual  middle  sensor  if 
m,  the  number  of  array  elements,  is  odd,  otherwise  it  is  a  fictitious  center.  When  m  is  even, 
the  choice  of  the  fictitious  center  as  the  coordinate  reference  corresponds  to  setting  p  =  y 
and  r  =  |  in  the  general  case. 

For  a  one-to-one  correspondence  between  the  values  of  <fr:  and  tvj,  we  shall  assume  that 
d  <  A/2  (reference  [18],  page  27).  Thus,  by  reorder  the  two  plane  waves  if  necessary, 
without  loss  of  generality,  we  can  always  write  the  direction  matrix  as 


A(origin  at  Q)  = 


e-i(p-l+r)u/ 

e-i(p-l+r)(u>+Au>) 

e-i(l+r)w 

e-i(l+r)(w+Aw) 

e-irw 

e-ir(w+Au<) 

ei(l -r)u< 

ei(l— r)(o>+Au>) 

ei  (m-p-r)w 

gi(m— p-r)(ui+Au/) 

(11) 
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where  w  and  w  +  Aw  are  electrical  phase  angles  with  Aw  <  ir.  Consequently, 

A(origin  at  Q)  =  e'i(P-1+’>  U  A0(origin  at  Q),  (12) 


where 


U  = 


-i2  u 


gi(m-l)w 


,  Ao(origin  at  Q)  = 


1  e-i(p-l+r)Aw 
1  e-*(l+r)Au» 

1  g— irAw 

1  ei(l-r)Aw 

j  gi(m— p— r)Au> 


•  (13) 


Note  that  U  is  unitary,  that  is  UU*  =  I  =  U*U. 
In  the  same  manner, 


D(origin  at  Q )  = 


_j(p  _  1  +  r)e-i(p-l+r)u-  _j (p  _  l  +  r)e-i(p-l+r)(u;+Au;) 


-i(l  +  r)e'i(1+r)‘‘ 


-ire 


i(l  -  r)ei<1-r)w 


-i(l  +  r)e-i(i+r)(^+Au,) 

_ire_ir(u'+Aw) 
i(l  -  r)ei(l-r)(u-+AW) 


can  be  factored  as 


where 


i(m  -  p  -  r)ei(m-P-’>  i(m  -  p  -  r)ei(m~P-r)("+Au») 


D(origin  at  Q)  =  e~^-^r)w  U  D0(origin  at  Q), 


-(p  -  1  +  r)i  -(p  -  1  +  r)ie-^-1+r)^ 


D0(origin  at  Q)  = 


— i(  1  +  r) 
-ir 

i(l  -r) 


-i(l  +  r)e-i<1+r)Au 
-i  re~irAw 
i(l  -  r)ei<1-r>A" 


(m-p-r) i  (m-p-  r)ie^m~P-r^ 


(14) 


(15) 


(16) 
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Using  expressions  (12)  to  (13)  and  (15)  to  (16)  together  with  the  fact  that  U  is  unitary, 
we  see  immediately 

(A*A)-1  =  (ASAo)'1,  (17) 

A+  =  A+IT,  (18) 

Pa  =  UPAoU%  (19) 

PA  =  UPA0U‘,  (20) 

H  =  DJP^Do,  (21) 

R  =  URqU*  (22) 

where  Ro  =  A0SA|5  +  a2I. 

The  expression  (22)  states  that  R  and  Ro  are  similar.  It  is  well  known  that  similar 
matrices  have  the  same  eigenvalues  with  the  same  multiplicities.  Also,  from  (22),  it  follows 
that  if  f  is  an  eigenvector  of  Ro  then  Uf  is  an  eigenvector  of  R  corresponding  to  the  same 
eigenvalue,  and  if  e  is  an  eigenvector  of  R  then  U*e  is  an  eigenvector  of  Ro  corresponding 
to  the  same  eigenvalue.  Thus, 

E,  =  UFS  (23) 


where 


with  the  columns  fi ,  •  •  • ,  fj  of  F,  are  normalized  eigenvectors  corresponding  to  eigenvalues 
Aj, . . . ,  Xd  of  Ro  (also  of  R).  Apply  (18),  (21)  and  (23)  to  (7)  and  (8)  we  have 

V"  =  -2Re{H  0  (Aq  FsWF*Aq  *)T}  (24) 

Q  =  2<72Re{H0(A+FsWAsA-zWF;A+*)r}.  (25) 

That  is,  in  the  case  of  two  narrowband  plane  signals  impinging  on  a  uniformly  spaced  linear 
array,  the  asymptotic  covariance  matrix  C wsf  of  the  estimation  error  using  subspace  fitting 
with  the  weighting  matrix  W  in  Theorem  2  is  dependent  on  Aw  and  independent  of  w. 
For  two  arrivals  with  signal  powers  tt2  >  0,  the  source  covariance  matrix  S  is  given 

by 


TTi  \/?rl7r2  P 

V^l  7T 2  P  7T2 


ys  w 

w-"  v^J’ 
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(26) 


where  p  is  the  normalized  correlation  between  the  first  and  second  baseband  signals  and 
0=arg(p)  is  the  temporal  phase  difference.  If  p= 0,  the  two  signals  are  incoherent,  \p\=\, 
they  are  coherent,  and  0  <  |p|  <  1,  they  are  partially  correlated.  Since  el6=  e'(5+2 *»)  for  any 
integer  k,  we  shall  consider  only  0  <  0  <  2x.  Note  that  0  depends  on,  as  being  measured 
with  respect  to,  the  coordinate  reference  of  the  array. 

Apply  (17)  and  (21)  to  (2)  to  (5)  and  we  see  readily  that  CRB  and  C mu  are  dependent 
on  A (jj  and  independent  of  u. 

We  summarize  the  above  discussion  in 

Lemma  1  For  any  two  distinct  narrowband  plane  signals  with  electrical  angles  u \  and 
u>2  impinging  on  a  uniform  linear  array,  for  any  choice  of  coordinate  origin  along  the  line 
segment  connecting  the  first  and  last  sensing  elements,  the  asymptotic  covariance  matrices 
CRB,  Cmu,  Cwsf  (thus,  in  particular,  Cml  )  are  dependent  on  Aw=|wi  -  W2I  arid  inde¬ 
pendent  of  U\  and  u>2-  Furthermore,  the  result  holds  for  arbitrary  signal  correlation  including 
full  coherence,  except  for  C mu- 

Since  e±‘(fl+,r)  =  -e±w,  by  applying  (26)  to  (2),  we  obtain  immediately  the  following 
result  which  can  be  deduced  from  references  [13],  [14]. 

Lemma  2  For  any  coordinate  reference  of  the  array  along  the  line  connecting  the  first  and 
last  sensors,  the  elements  of  the  2x2  real  symmetric  matrix  CRB,  considered  as  functions 
of  the  temporal  phase  difference  0,  satisfy 


[CRB(0 +  *)],, 

=  [CRB(0)]t<  for  i  =  l,2 

(27) 

[CRB(0  +  tt)]12 

=  -[CRB(0)]12. 

(28) 

The  matrix  H  is  a  prominent  component  of  CRB,  Cmu ,  and  Cwsf  as  seen  in  (2),  (4), 
(7)  and  (8),  as  well  as  in  (9).  The  next  result,  whose  proof  appears  in  Appendix  A,  gives  a 
general  description  of  H. 
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Lemma  3  For  any  choice  of  coordinate  origin  on  the  line  segment  connecting  the  first  and 
last  array  elements,  H  is  of  the  form 


H  = 


|H12|  <  Hn,  Hu  >  0. 


Hu  H12 
h12  H„ 

Furthermore,  let  Q  be  the  general  coordinate  reference  given  in  (10)  then 


Hu  (origin  at  0) 

m(m  -  l)(2m  —  1) 


rn2  - 


m  j 

/sin (a^)\2  j 

V  J 


-  m(p  -  1  +  r)(m  -  p  -  r) 
m2  (m  +  1  -  2 (p  4-  r))2 


-  (m  +  1  -  2(p  +  r))  Re  ^fc(origin  at  Q)  5Z(origin  at  O)) 

+|Efe(°riginatO)|2}, 


and 


Hi2(origin  at  O) 

=  (origin  at  O) 
1 


|m2  (m  +  1  -  2(p  +  r))  ^Tjorigin  at  Q) 


a  f“lf>y 

V  ««(¥) ) 

^(m_+l  -2(p  +  r)f  E(originit0) 

-  (E((ciriSin  at  O))2  H(°riSin  at  0)}  • 


where 


m— 1 


(29) 


(30) 


^(origin  at  Q)  ^  E  «“ =  (^iy)  (31) 


Et(origiaatO)=  E(*-(J>-1  +  ’-))ei(M’',+r,‘1"  =  -i^(E('>'iSinatO)),  (32) 


Jk=0 

m-1 


53^ (origin  at  O)  =  £(fc  -  (p  -  1  +  r))2ei(*  (p  1+r))Aw  =  -i—  (£fc(origin  at  O))  • 

(33) 


11 


For  ease  of  notation  we  define 


A 

a  = 


m 


2  _ 


m 


/■n(afttt) 

V 


) 


2’ 


where 


Notice  that 


p(Aw;m)  = 


sin(sf*) 

«n(^)  ' 


m 


2 


( 


sin(^)  J 


2 


=  det(A*A)  =  det(ASAo) 


>  0 


(34) 


(35) 


(36) 


since  AJJAo  is  Hermitian  (hence  its  determinant  is  non  negative)  and,  for  0  <  Aw  <  tt  and 
m  >  3,  the  two  columns  of  Aq  are  linearly  independent  (hence  det  (AJAq)  ^  0).  Thus, 


sin(gfa) 

sin(^) 


<  m. 


Therefore, 


a  >  |6|  >  0,  sign  (6)  =  sign  (p  (Aw;  m)) . 


(37) 


As  a  consequence  of  Lemma  3,  if  the  coordinate  origin  is  at  the  center  of  the  array 
(hence  p  +  r  =  regardless  of  the  parity  of  m),  it  follows  from  (31)  to  (33)  that  H  is  a 
real  matrix  with  entries 

,,  .  .  .  .  m(m  -  l)(m  +  1)  (  dp  \2 

Hn  (origin  at  center)  =  — - - — - “  VdAwj  ’  (38) 

and 

H.aforigin  at  center)  =  -  jAiL;  -  b  ■  (££)  .  (39) 

Evidenced  by  (2),  (4),  (7),  and  (8),  as  well  as  in  (9),  the  choice  of  the  center  of  the 
array  as  the  coordinate  reference  will  greatly  facilitate  the  investigation  of  the  effects  of  the 
temporal  phase  difference  on  CRB,  C mu,  and  C wsf  since  the  cross  terms  that  involve 
both  the  temporal  phase  difference  9  and  arg(Hn)  (dependent  on  the  electrical  angular 
separation  Aw)  will  be  simplified  when  H  is  real.  Clearly,  the  analysis  can  be  carried  out  in 
terms  of  the  general  coordinate  reference  (10).  However,  to  avoid  cumbersome  arguments 
and  expressions  involving  awkward  translations,  we  elect  to  present  results  in  the  remainder 
of  the  paper  for  the  coordinate  origin  chosen  at  the  center  of  the  array. 
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4.  EFFECT  OF  TEMPORAL  PHASE  DIFFERENCE  ON  CRB 


The  main  result  in  this  section  can  be  summarized  as  follows.  Recall  that,  by  Lemma 
2,  we  can  restrict  the  temporal  phase  difference  9  to  the  interval  [0,  ?r] . 

Proposition  A  Consider  [CRB]it ,  i  =  1,2,  as  a  function  of  the  temporal  phase  dif¬ 
ference  9,  0  <  9  <  n,  of  two  narrowband  signals  impinging  on  a  ULA  whose  coordinate 
origin  is  chosen  at  the  center  of  the  array.  Then  [CRB(0)]„  ,  i  =  1,2,  is  symmetric  about 
tv/2.  Furthermore,  it  is  decreasing  for  0  <  9  <  t/2,  increasing  for  7r/2  <  9  <  ir,  attains 
the  maximum  (either  finite  or  infinite)  at  0  and  n,  and  the  minimum  at  tt/2. 


Detailed  statements  and  proofs  leading  to  this  result  are  presented  below  in  a  sequence 
of  two  lemmas.  Proposition  A  can  be  deduced  from  reference  [15];  however,  our  argument 
is  based  on  expression  (2)  derived  in  reference  [2]. 

First,  using  the  coordinate  origin  at  the  center  of  the  array,  we  obtain  from  Lemma  3 


and  then 


Re(H  0  ST)  =  y/iriTti 


yj¥2  Hu  \p\  H12  Re(e~'6) 
\p\  H12  Re(ew)  ^Hu 


det  (Re(H  0  ST))  =  ttjttj  (hu2  -  |p|2  Hi22  cos2(0))  . 


(40) 


We  shall  consider  [CRB],-,- ,  i  =  1,2,  as  a  function  of  9.  In  view  of  Lemma  2,  it  suffices 
to  restrict  9  to  the  interval  0  <  9  <  ir. 

Using  (2)  and  (40),  we  have 


[CRB]ti 


Clearly, 


*'=1,2. 


(41) 


Lemma  4  Consider  [CRB]tt ,  i  =  1,2,  of  two  narrowband  plane  waves,  with  electrical 
angular  separation  Au>,  impinging  on  a  ULA  whose  coordinate  origin  is  chosen  at  the  center 
of  the  array.  Suppose  that  either  the  two  signals  are  incoherent  (p  =  0)  or  Hi2  =  0,  then 

2 

[CRB]»  =  2^l7’  *  =  1’2'  <42) 
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that  is  the  variances  [CRB]n  and  [CRB]22  are  independent  of  9. 

We  note  that,  using  (2)  and  (26),  Lemma  4  holds  for  arbitrary  coordinate  reference 
where  Hu,  Hi2  and  [CRB]-,  i  =  1,2,  are  measured  with  respect  to  that  coordinate 
reference.  Also,  using  (30)  to  (33),  we  observe  that  for  any  choice  of  coordinate  reference 
and  for  any  even  m  >  4,  then  H12  =  0  when  Aw  =  y  BW  (=  ir)  where  BW  =  ^  denotes 
a  beamwidth. 

Next,  we  have 


Lemma  5  Consider  [CRB],-,-,  *  =  1,2,  as  a  function  of  9,  for  0  <  9  <  n.  Then 
[CRB(0)]„  ,*  =  1,2,  is  symmetric  about  | .  Moreover ,  [CRB(0)],-,- ,  i  =  1,2,  is  decreasing 
for  0  <9  <  x/2,  increasing  for  x/2  <  9  <  tt,  attains  the  maximum  at  0  =  0,  7 r,  and  the 
minimum  at  9  —  ir / 2,  with 

min[CRB(«)]„  =  ^,  <  =  1,2,  (43) 


max[CRB(0)]„ 


2xi  Hu 


i  =  1,2. 


Furthermore,  [CRB],,- ,  i  =  1,2,  is  arbitrarily  large  (00)  if  and  only  if 


(44) 


(i)  |Hi2|  =  Hu  (hence,  H  is  singular), 

(ii)  |p|  =  1  (hence,  S  is  singular  i.e.,  the  two  signals  are  coherent),  and 

(iii)  9  —  0,  7r  (i.e.,  the  two  signals  are  in  phase  or  out  of  phase). 


Proof.  In  view  of  Lemma  4,  without  loss  of  generality,  we  can  assume  0  <  |^|  <  1  and 
H12  ^  0.  We  notice  that  for  any  £  €  [0, 7r/2], 

c°s2  (|  ±0)  =  sin2(0- 

Hence,  for  0  <  9  <  tt,  the  variance  [CRB(0)]„ ,  t  =  1,2,  is  symmetric  about  y.  The  rest 
of  the  lemma  follows  from  examining  the  factor 


1 


1 
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of  [CRB],-,- ,  i  =  1,2,  in  (41)  as  a  function  of  9.  M 


Qualitatively,  when  the  normalized  correlation  magnitude  |/>j  is  large  (close  to  unity) 
then  it  follows  from  (41)  that,  other  things  being  equal,  effect  of  9  on  [CRB]^  ,  i  =  1,2,  is 
intensified. 

On  the  other  hand,  when  the  electrical  angular  separation  Au;  is  small  in  terms  of 
beamwidth,  (equivalently,  the  two  steering  vectors  which  are  columns  of  A  are  close),  then 
the  two  columns  of  D  (being  elementwise  derivatives  of  the  two  steering  vectors)  are  close. 
Since 

H  i  D'PiD  =  D-(Pi)‘pjjD  =  (PiD)'PiD, 

it  follows  that  the  ratio  |Hi2|/Hn  is  close  to  unity  when  Au;  is  small.  Also,  from  (38)  we 
notice  that  Hu  is  small  when  A u  is  small.  Consequently,  other  things  being  equal,  from 
(41)  we  see  immediately  that  effect  of  9  on  [CRB],-,- ,  i  =  1,2,  is  augmented  when  Au;  is 
small. 

Therefore,  the  effect  of  temporal  phase  difference  on  [CRB]ti ,  i  =  1,2,  are  intensified 
when  the  two  signals  are  closely  spaced  and/or  their  normalized  correlation  magnitude  is 
high. 

For  the  general  coordinate  reference,  we  notice  that  H  need  not  be  real  and  thus  (41) 
is  replaced  by 


[CRB],,  = 


27r«H“  \l-jp|2  (lgiil)  cos2(arg(H u)  -  9) 


,  i  =  1,2.  (45) 


In  Appendix  B,  we  present  an  example  where  the  electrical  angular  separation  Aw  is  an 
integer  multiple  of  (standard  or  first  null)  beamwidth,  denoted  by  BW,  with  one  BW  equal 


to  2n/m.  That  is 

Au;  =  l  BW  =  l  , 

where  l  =  1,2, ...,  [yj  with  [yj  denotes  the  largest  integer  less  than  or  equal  to  y.  It 
is  interesting  to  note  that,  in  this  situation,  H  is  nonsingular  for  m  >  4  and  singular  for 
m  =  3.  Therefore,  for  m  >  4  and  Au;  =  l  BW,  t  =  1,2,...,  [yj,  the  variance  [CRB],,  , 
i  =  1,2,  is  finite  for  arbitrary  signal  correlation,  whereas,  when  m  =  3  and  Au;  =  1  BW, 
the  variance  [CRB],, ,  i  =  1,2,  is  arbitrarily  large  when  the  two  signals  are  correlated  with 
an  appropriate  temporal  phase  difference  9  depending  on  the  coordinate  reference. 
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5.  EFFECT  OF  TEMPORAL  PHASE  DIFFERENCE  ON  CMU 


With  the  coordinate  reference  at  the  center  of  the  array,  by  the  computation  in  Appendix 
C  we  have 

1CmuI"  =  2^  ii77(i  -  IpF)  ' +  Bi  cos{>))  (46) 

where 


M  =  1  +  (*,  +  \p\2  it,)  •  a,  (47) 

Bi  =  2v/7ri*-2  jp| .  b,  (48) 


with 


Hn  given  in  (38), 
a  and  6  given  in  (34), 

det(S)  =  7T17T2  (l  -  |p|2)  >  0  and  i,j  £  {1,2},  j  ^  t. 


Notice  that  A{  and  Bi,  i  =  1,2,  are  independent  of  the  temporal  phase  difference  6. 
Also,  using  the  fact  that  2y'jr1x2  \p\  <  *j  +  \p\2  *r«,  t,j  =  1,2,  j  i  (as  “geometric  mean” 
<  “arithmetic  mean”  or,  simply,  by  a  quadratic  expansion)  in  conjunction  with  (37),  we 
obtain,  for  t  =  1,2, 


A,  >  |£,|  >  0,  sign  (Bi)  =  sign(p(Aw;m)) . 
If  p  =  0  then  Bi  =  0  and,  using  Lemma  4,  we  obtain  immediately 


Proposition  B.l  Consider  [Cjvfc/]tj »  »  =  1,2,  of  two  narrowband  plane  waves,  with 
electrical  angular  separation  Aw,  impinging  on  a  ULA  whose  coordinate  origin  is  chosen 
at  the  center  of  the  array.  Suppose  that  the  two  signals  are  incoherent  (p  =  0)  then,  for 
*  =  1,2, 


2»iHu  '  (*  +  “ 
mjn[CRB(0)]„  +  2  Hu  •  a  •  (min  [CRB(0)]*  j  . 
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Next,  we  consider  the  general  case  (0  <  \p\  <  1)  where  Proposition  B.l  depicts  a 
degenerate  situation.  For  any  0  <  Au  <  ir,  there  exists  a  nonnegative  integer  l,  with 
0  <  i  <  [yj,  such  that  £BW  <  Au>  <  min{(£  +  l)BW,7r};  where  BW  =  2w /m  denotes  a 
beamwidth,  and  [xj  denotes  the  largest  integer  less  than  or  equal  to  x.  Clearly,  for  *  =  1,2, 

Bi  >  0,  for  £BW  <Aw<(£  +  l)BW  and  £  even  (including  zero), 

B{  <  0,  for  l  BW  <  Au;  <  (£  +  1)  BW  and  £odd, 

Bi  =  0,  for  Au  =  (£  +  1)  BW. 

Therefore, 


Proposition  B.2  Consider  [Cm£/]„  ,  i  =  1,2,  as  a  function  of  the  temporal  phase 
difference  9,  0  <  9  <  2n,  of  two  narrowband,  noncoherent  (0  <  |p|  <  1)  plane  waves,  with 
electrical  angular  separation  Alj,  impinging  on  a  ULA  whose  coordinate  origin  is  chosen  at 
the  center  of  the  array.  Suppose  t  BW  <  Au  <  nain{(£  +  1)  BW,7r}  for  some  nonnegative 
integer  l  with  0  <  l  <  [yj.  Then  [CmlK#)];, ,  *  =  1,2,  is  symmetric  about  9  =  ir. 
Moreover, 

(a)  for  f  BW  <  Au  <  min{(£  +  1)  BW,  jt}, 

•  if  l  is  even  (including  zero),  then  [Cwy^))-  ,  i  =  1,2,  is  decreasing  for 
0  <  9  <  7r,  increasing  for  tt  <  9  <  27T,  attains  its  maximum  at  9  =  0,  2tt, 
and  its  minimum  at  9  =  ir-, 

•  if  l  is  odd,  then  [Ca/1/(0)]„  ,  *  =  1,2,  is  increasing  for  0  <  6  <  ir, 
decreasing  for  n  <  9  <  2ir,  attains  its  minimum  at  9  =  0,  2ir,  and  its 
maximum  at  9  =  ir, 

in  either  case, 

-rtCiwWk  =  2TjHlT(i-|,|2)'<i4'  '  i  =  1’2 

=  2TiH,T(r-WJ)'('4i  +  lal)’  i=1,2 

where  Ai  and  Bi  are  given  in  (47)  and  (48); 


(b)  for  Au>  =  (£+1)BW(<tt), 


[Cmc/W],,  = 


2 *»  (1  -  H2)  HiiLh>=(*+i)bw 


■|4»Iau»= 


(<+l)BW 


*  =  1,2, 
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where 


^■Iaw=(<+i)bw  1  +  mdet(S) 


(»j  +  \P?  *.)  i  »,j  =  1,2  and  j  ±  i. 


and 


H, 


m 


llllAw=(/+l)BW  - 
hence,  it  is  independent  of  0. 


m2  -  1  1 


3  sin’(^)J 


Furthermore,  in  all  cases,  when  \p\  approaches  unity,  [Cm[/(#)];,  gets  arbitrarily  large, 
regardless  of  the  temporal  phase  difference  6  and  electrical  angular  separation  A w  or  any 
other  parameters. 


It  follows  immediately  from  (46)  to  (48)  that  for  large  |p|,  other  things  being  equal, 
effect  of  0  on  (CMt/]j; ,  t  =  1,2,  is  intensified. 

On  the  other  hand,  when  Aw  is  small  (in  terms  of  beamwidth),  Hu  is  small  and,  by 
(34),  a  and  6  are  large.  Therefore,  effect  of  0  on  [CAfi/],j ,  t  =  1,2,  is  intensified  when  Aw 
is  small. 

Consequently,  effect  of  the  temporal  phase  difference  on  [Cjvft/],, ,  i  =  1,2,  is  intensified 
when  the  two  signals  are  closely  spaced  and/or  their  normalized  correlation  magnitude  is 
high. 
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6.  EFFECT  OF  TEMPORAL  PHASE  DIFFERENCE  ON  CML 


From  (9),  with  the  coordinate  reference  at  the  center  of  the  array,  a  straightforward, 
but  tedious,  computation  yields 


\CML)ti  =  [CRB],,  +  2  H„  •  a  •  ([CRB],,2  +  [CRB]122) 

-  4  H,2  •  b  ■  [CRB]ti  •  [CRB]12  ,  i  =  1,2,  (49) 


where 


[CRB],, 

[CRB]12 


<T2 

r  1  ' 

27 r,  Hn 

-<7*\p\H 

l1-! 
12  1 

pi2  (is) 

f 

cos2(0) j 

:os (8) 

> 

\ 

2v/7r17T2  H 

2 ' 

11  1 

{i  -IpI2  ( 

m ) 

(50) 

(51) 


Clearly,  if  either  p  =  0  or  H12  =  0,  then  [CRB(0)]12  =  0  and  [CRB(0)]tl  = 
min#  [CRB(0)J- ,  i  =  1,2.  Therefore, 


Proposition  C.l  Consider  [Cjutli,  >  1  =  1.2,  of  two  narrowband  plane  waves,  with 
electrical  angular  separation  Aw,  impinging  on  a  ULA  whose  coordinate  origin  is  chosen  at 
the  center  of  the  array.  Suppose  that  either  the  two  signals  are  incoherent  ( p  =  0)  or  Hi2 
=  0  then,  for  i  =  1,2, 

[C ml]«  =  min[CRB(d)],  +  2  H„  •  a  •  (mjn  [CRB(^)],,)2  , 
which  is  independent  of  8  and  equal  to  [Cjvfr/],,  when  p  =  0. 


Proposition  C.l  is  consistent  with  the  general  result  of  Stoica  and  Nehorai  (reference 
[2]),  which  states  that  if  the  matrix  S  (not  necessarily  a  2  x  2  matrix)  is  diagonal  then  the 
asymptotic  covariance  of  ML  coincides  with  that  of  MUSIC. 

We  next  consider  the  general  case  (0  <  \p\  <  1  and  0  <  |Hn|  <  Hu)  where 
Proposition  C.l.  is  a  degenerate  situation. 

First,  we  suppose  that  0  <  |pj  <  1.  To  facilitate  the  discussion,  we  rewrite  the  expression 
(49)  as  follows. 

[CM  Mi  =  Tx{8)  +  T2(8)  (52) 


19 


where 


Ti(9) 

T2(9) 


=  [CRB(0)].,  +  2  H „  •  a  •  ([CRB(fl)]f ,2  +  [CRB(0)]122)  , 


<r4  HH122 


cos(0) 


^  \r\  ^  _ 

'(l-|/.P(fe)2cosW) 


2  • 


(53) 

(54) 


Notice  that 


dTt 

dO 


-  sin(0) 


7T,  Ha 


(i-W  (fe)2“*2w)' 


x  F{6)  cos (8), 


(55) 


Obviously,  if  either  p  =  0  or  Hi2  =  0  then  both  F{9)  and  G{8 )  vanish,  and  hence 
[Cmz,(0)]„  *s  independent  of  8  as  known  in  Proposition  C.l. 

For  p  ^  0  and  HJ2  /  Owe  have  immediately 


F(8)  >  0,  sign(G(0))  =  sign  (p(Au>;m)) . 


Thus,  T\{9)  is  positive  and,  as  a  function  of  9,  behaves  exactly  like  [CRB(0)]„  .  The 
term  T2{9)  of  (52)  is  more  troublesome  due  to  the  factor  p{Au\m). 

Recall  that,  for  any  0  <  Au  <  ir,  there  exists  a  nonnegative  integer  t,  with  0  <  l  <  [yj, 
such  that  £BW  <  Au>  <  min{(f  +  l)BW,7r};  where  BW  =  2 n/m  denotes  a  beamwidth 
and  [xj  denotes  the  largest  integer  less  than  or  equal  to  x.  Clearly, 


p(Auqm) 


>  0,  for  f  BW  <  Au  <  ((  -f  1)  BW  and  l  even  (including  zero), 
•  <  0,  for  ^BW  <  Au  <  (1+  1)BW  and  ^odd, 

=  0,  for  Au  =  (*+  1)BW. 
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Suppose  that  p(A u;m)  >  0.  Clearly,  the  global  maximum  occurs  at  9  =  0,  2x.  From 
(55)  and  (56)  we  have 


thus 


dTi 

dd 


dTi 

de 


<  0, 

0  <  9  <  tt/2,  and  x  <  9  <  3x/2, 

=  0, 

9  =  x/2,  x,  3x/2,  2x, 

>  0, 

x/2  <  9  <  x,  and  3x/2  <  9  <  2 x, 

<  0, 

0  <  9  <  x, 

=  0, 

0  =  0,  x,  2x, 

>  0, 

x  <  0  <  2x, 

^  [CWl(0)]« 


dTi  dTi 

de  +  do 


f  <  0,  0  <  9  <  tt/2, 

|  >  0,  3x/2  <9  <2i r, 


which  implies  [C ml(^)]^  is  decreasing  for  0  <  0  <  x/2  and  increasing  for  37r/2  <  6  <  27r. 

Now  consider  x/2  <  9  <  3x/2.  By  the  symmetry  of  [Cml(0)]„  about  7r  we  need  only 
examine  the  case  where  tt/2  <  9  <  x. 


We  write 


nHn 


X  (F(9)  cos (9)  +  G{9)) . 


Suppose  that  jg  ([0^(9)],,-)  is  not  identically  zero  on  (7r/2, 7r] ,  that  is  p  ^  0  and  Hu  ^  0. 
If  j^([Cj ml(9)\h)  vanishes  at  any  point  9q  E  (r/2,  x),  then  9q  is  a  zero  of  the  function 
F{9)  cos (0)  -f  G(9)  and 

IGWl  J_  n™  <  , 

F(9)  |p|  denom,  - 

where 


and 


denom  =  X1X2  ^1  -  |p|2  cos 2(9)'j 

+  ct2  •  a  •  ^2?rj  +  i r,  ^1  +  |p|2  cos2(9)^jJ  . 
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Notice  that  for  either  small  \p\,  or  small  |Hi2|/Hn  or  for  0  being  close  to  7r/2,  then 


num  A  o2  »/7ri7r2  •  |&|  1 

-  'y  —  - — -  <  — 

denom  *1*2  +  <r2  ■  a  •  (2  nj  4-  tt;)  2 


(57) 


sin(*f*) 

sin(^) 


<  m, 


7Ti  +  7T2 


Therefore,  if  ^  is  close  to  zero  then  jg  ([Cml(^)],,)  vanishes  at  some  6q  6  (ir/2,7r) 
which  is  close  to  tt/2.  The  case  where  the  normalized  correlation  magnitude  \p\  is  large  and 
the  SNR’s  ^  and  are  high  belongs  to  this  situation. 

On  the  other  hand,  0  <  \p\  <  7  then  $g  ([Cmz,(0)D  does  not  vanish  on  (zr/2, 7r) .  Since 


-^([CmlWU  <  0,  ^(|CMt(«)]„)  =  0 

dv  e=n/2  dd  0=* 

and  since  jg  ([CmlW];;)  is  a  continuous  function  in  6 ,  it  follows  that  jg  ([0^(0)]^)  <  0 

on  (7r/2, 7r) .  Consequently,  [CmiX#)]-  is  decreasing  on  (tt/2, zt)  and,  by  its  symmetry  about 

it,  increasing  on  (n,3n/2)  with  the  minimum  achieved  at  8  =  7r. 

Thus,  the  behavior  of  [Ca/£,(0)]^  on  (7r/2,  37t/2)  is  quite  intricate  due  to  the  interaction 

of  different  parameters.  In  general,  since 


^([C«my  <  0,  »e  (0,1/2], 

3j([c»«,(«)y|^  =  0, 

and 

^([CuMi)  >  0,  0e[37T/2,27T), 

and  since  ([CmiX#)],;)  attains  its  (global)  maximum  at  9  =  0,2x,  it  follows  that  there  is  a 
(global)  minimum  at  do  (E  (7r /2,  tt]  and,  by  its  symmetry  about  tt,  also  at  2t  —  90.  The  exact 
location  of  80  however  depends  of  Aw  and  other  parameters  tti/ct2,  |p|,  and  m. 

A  similar  argument  can  also  be  applied  to  the  case  where  p(Aw;m)  <  0. 

Finally,  we  note  from  (49)  (also  from  (9))  that,  for  0  <  \p\  <  1,  the  variance  [Ca/l],, , 
i  =  l,2,  can  be  expressed  in  terms  of  elements  of  the  matrix  CRB  which  hold  for  0  <  IpI  <  1. 
Therefore,  the  following  result,  which  is  a  summary  of  the  preceding  discussion,  also  holds 
when  |p|  =  1  by  an  elementary  continuity  argument. 


Proposition  C.2  Consider  [Ca/l];,- »  i  =  1,2,  as  a  function  of  the  temporal  phase  dif¬ 
ference  6,  0  <  9  <  2tt,  of  two  narrowband  plane  waves  with  arbitrary  signal  corelation 


(0  <  |p|  <  1)  and  electrical  angular  separation  Aw,  impinging  on  a  ULA  whose  coordinate 
origin  is  chosen  at  the  center  of  the  array.  Suppose  £  BW  <  Aw  <  min{(£  +  l)BW,7r} 
for  some  nonnegative  integer  £  with  0  <  l  <  [yj .  Then 

(a)  for  £  BW  <  Aw  <  min{(i  +  l)BW,5r},  then  [Cml (0)]„  »  *  =  1,2,  is  symmetric 
about  6  —  it  and 

•  if  £  is  even  ( including  zero)  then,  on  the  interval  [0,  *r],  [Cml(Q)\u  attains 
its  (global)  maximum  at  6  =  0 ,  and  its  (global)  minimum  at  6o  £  (*r/2,  jr]; 

•  if  £  is  odd  then,  on  the  interval  [0, 7r],  [0^(0)],, ,  i  =  1,2,  attains  its 
(global)  maximum  at  6  =  ir,  and  its  (global)  minimum  at  Oq  €  [0, 7r/2); 

In  either  case,  the  exact  location  of  0o  depends  on  Aw  and  other  parameters  iti/cr2,  it 2/ff2, 
|p|,  and  m; 


Furthermore,  in  all  cases,  ,  t  =  1,2,  is  arbitrarily  large  (00)  if  and  only  if 
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(i)  |Hi2 1  =  Hu  (hence,  H  is  singular), 


(ii)  \p\  =  1  (hence,  S  is  singular  i.e.^the  two  signals  are  coherent),  and 

(iii)  9  =  0,  ir  (i.e^the  two  signals  are  in  phase  or  out  of  phase). 

Qualitatively,  it  follows  from  (49),  using  (50)  and  (51),  that  effect  of  the  temporal  phase 
difference  on  ,  i  =  1,2,  is  intensified  when  the  normalized  correlation  magnitude  jpj 

is  high. 

On  the  other  hand,  since 


,  i  =  l,2, 


H„  •  [CRB]'  =  -4—  - — 

V  1  -  I/>12  (S 


\p?  (h^)  cos2W 


HX1  [CRB&  = 


\p\2  /h« 

47TXX2HH  \Hu 


cos  (9) 

VHni  (l-l/Wfe)2 


‘  4tjT2  Hu  VH„7  ’ 

•  [CRB],,  ■  [CRB],,  =  •  (|^)2  ■  - - 

'^H"  vh„;  (j _  w  («*)’  COi 


and  since  Hu  is  small  and  a,  b  are  large  when  A u>  is  small  (in  terms  of  beamwidth),  it  follows 
from  (49)  that,  other  things  being  equal,  effect  of  9  on  [Ca/l]„  ,  t  =  1, 2,  is  augmented  when 
the  two  signals  are  closely  spaced. 

In  short,  effects  of  the  temporal  phase  difference  on  [Ca/l],;,  i  =  1,2,  are  intensified 
when  the  two  signals  are  closely  spaced  and/or  their  normalized  correlation  magnitude  is 
high. 
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7.  NUMERICAL  RESULTS 


For  simplicity,  we  shall  consider  the  situation  of  two  plane  waves  of  equal  powers  im¬ 
pinging  on  a  uniform  linear  array.  By  setting  ni  =  W2  in  (41),  (46)  to  (48)  and  (49)  to  (51), 
we  have 

[CRBjjj  —  [CRB]22  i  =  [C\fl/]22  *  =  [^Aft] 22  •  (^8) 

This,  however,  can  be  obtained  without  actually  computing  the  corresponding  variances  as 
in  Appendix  A. 

In  figures  1(a)  to  3(d),  we  consider  (square-roots  of)  [CRB],,-,  ,  and  [Cml],; 

as  functions  of  two  variables,  the  electrical  angular  separation  Au>  and  the  temporal  phase 
difference  6  where 

0.1  BW  <  Aw  <  1  BW,  and  0°  <  9  <  360°. 

The  vertical  axis  of  each  three-dimensional  figure  depicts  y/N  times  of  the  corresponding 
standard  deviation  (see  (3)),  measured  in  BWs,  where  N,  the  number  of  time  samples 
(snapshots),  is  assumed  to  be  large. 

Furthermore,  all  figures  have  the  following  common  parameters 


m  =  5, 


7T2 


=  20  dB. 


The  first  set  of  figures  (figures  1(a)  to  1(d))  presents  ^/[CRB]--,  for  \p\  =  0.50,  0.90, 
0.95,  and  1.00,  respectively.  The  second  set  (figures  2(a)  to  2(d))  presents  \J{Cmu\h  for 
|p|  =  0.50,  0.90,  0.95,  and  0.97,  respectively.  Finally,  is  presented  in  the  third 

set  (figures  3(a)  to  3(d))  for  |p|  =  0.50,  0.90,  0.95,  and  1.00,  respectively. 

Figures  1(a)  to  1(c),  2(a)  to  2(d),  and  3(a)  to  3(c)  are  plotted  using  general  formulas 
(2),  (4),  and  (9),  respectively;  figure  1(d)  is  based  on  (2)  and  figure  3(d)  on  (6)  together 
with  (7)  and  (8).  Thus,  this  graphical  sample  can  be  used  to  cross-examine  the  consistency 
of  our  results  with  respect  to  general  formulas  given  in  references  [2],  [3],  and  [6]. 
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Figure  1(d).  CRB  for  two  equipowered  signals 
with  |rho|-1.00,  SNR-20  dB,  impinging  on  a  5-element  ULA 
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Figure  2(c).  MUSIC  for  two  equipowered  signals 
with  |rho|=0.95,  SNR =20  <18,  impinging  on  a  5-element  ULA 


Figure  2(d).  MUSIC  for  two  equipowered  signals 
with  |rho|  =0.97,  SNR =20  dB,  impinging  on  a  5-element  ULA 
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Figure  3(a).  ML  for  two  equipowered  signals 
with  |  rho  1=0.50,  SNR =20  dB,  impinging  on  a  5-element  ULA 


Figure  3(b).  ML  for  two  equipowered  signals 
with  | rho  |  =0.90,  SNR-20  dB,  impinging  on  a  5-6iement  ULA 
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Figure  3(c).  ML  for  two  equipowered  signals 
with  |rho|  =0.95,  SNR=20  dB,  impinging  on  a  5-eiement  ULA 


Figure  3(d).  ML  for  two  equipowered  signals 
with  |rho| -1.00,  SNR-20  dB,  impinging  on  a  ULA 
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Notice  that  the  behavior  of  [CRB]|t ,  [Cjvft/],-,-  *  and  with  respect  to  the 

temporal  phase  difference  0  agrees  with  our  results  (Lemma  5,  and  Propositions  B.2  and 
C.2).  In  particular,  when  A w  =  1BW,  [C mu)u  is  independent  of  8  in  figures  2(a)  to  2(d) 
as  stated  in  Proposition  B.2,  and  also  when  Aw  =  1BW,  [Ca/l],-,-  behaves  in  the  same 
manner  as  [CRB]it  in  figures  3(a)  to  3(d)  as  expected  from  Proposition  C.2. 

Qualitatively,  as  discussed  in  sections  4  to  6,  effects  of  temporal  phase  difference  8 
on  [CRB]^  ,  [Cjv/(/]jj  and  [Ca*x] are  intensified  when  the  two  signals  are  closely  spaced 
and/or  when  their  normalized  correlation  magnitude  |p|  is  high.  Figures  1(a)  to  3(d)  exhibit 
this  behavior. 


8.  CONCLUSIONS 

We  derive  closed  form  expressions  for  the  Cramer-Rao  Bound,  MUSIC,  and  ML  asymp¬ 
totic  variances  corresponding  to  the  two-source  direction-of-arrival  estimation  where  sources 
are  modeled  as  deterministic  signals  impinging  on  a  uniform  linear  array. 

The  choice  of  the  center  of  the  array  as  the  coordinate  reference  results  in  compact 
expressions  that  greatly  facilitate  our  study  of  effects  of  temporal  phase  difference  (correla¬ 
tion  phase)  of  the  two  sources  on  asymptotic  variances  of  estimation  errors.  For  instance, 
the  behavior  of  the  Cramer-Rao  Bound,  which  was  obtained  earlier  by  Evans,  et  al.  (refer¬ 
ence  [15]),  can  now  be  stated,  with  respect  to  the  coordinate  reference  at  the  center  of  the 
array,  as  follows. 

(•)  The  CRB  variance  of  estimation  error,  considered  as  a  function  of  the 
phase  difference  8,  is  periodic  of  period  180°,  symmetric  about  90°,  and  on  the 
interval  where  0°  <  8  <  180°,  it  is  decreasing  when  0°  <  8  <  90°,  increasing 
when  90°  <  8  <  180°,  and  it  assumes  the  minimum  at  90°  and  maximum  (either 
finite  or  infinite)  at  0°  and  180°  (that  is,  either  the  two  signals  are  in  phase  or 
out  of  phase). 

Since  these  effects  are  intensified,  as  reflected  in  derived  expressions  and  exhibited  in 
numerical  results,  when  the  two  signals  are  closely  spaced  and/or  when  their  normalized 
correlation  magnitude  is  high,  the  following  is  a  summary  of  our  analytical  results  for 
MUSIC  and  ML  restricted  to  the  case  when  the  electrical  angular  separation  is  within  one 
beamwidth  (0  <  Aw  <  1  BW). 
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(•)  Suppose  0  <  Au;  <  1  BW,  then  the  asymptotic  MUSIC  variance  of 
estimation  error,  considered  as  a  function  of  the  temporal  phase  difference  9 , 
is  symmetric  about  180°,  decreasing  when  0°  <  9  <  180°,  increasing  when 
180°  <  9  <  360°,  attains  its  maximum  at  0  =  0°  (that  is,  the  two  signals  are  in 
phase)  and  minimum  at  9  =  180°  (that  is,  the  two  signals  are  out  of  phase). 

(•)  Suppose  0  <  Aw  <  1  BW,  then  the  asymptotic  ML  variance  of  estimation 
error,  considered  as  a  function  of  the  temporal  phase  difference  9 ,  is  symmetric 
about  180°,  and  on  the  interval  where  0°  <  9  <  180°,  it  attains  the  maximum 
at  9  =  0°  (that  is,  the  two  signals  are  in  phase)  and  minimum  at  9  =  9q  with 
0°  <  9  <  9o ,  where  the  exact  location  of  9q  depends  on  Au>  and  other  parameters 
(number  of  array  elements  m,  normalized  correlation  magnitude  |p|,  and  signal- 
to-noise  ratios  %  and  ^). 

Our  results  provide  an  analytical  background  for  simulations  such  as  those  performed 
in  references  [10]  and  [12]. 
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APPENDIX  A 


Proof  of  Lemma  3 

Although  the  main  purpose  of  this  appendix  is  to  prove  Lemma  3,  the  fundamental 
development  can  be  used  to  show  that  if  the  two  plane  waves  have  equal  powers,  that  is 
7T x  ==  fl"2,  then 

[CRB]n  =  [CRB]22  ,  [Cmi/]h  =  [Cmi/]22  ’  [cml]h  =  [Cml]22  >  (A-1) 

without  actually  computing  the  corresponding  variances. 

To  facilitate  the  discussion  we  introduce  the  following  family  of  2  x  2  Hermitian  matrices. 
Let 


that  is,  C  is  the  collection  of  2  x  2  Hermitian  matrices  M  with  Mu  =  M22-  The  following 
properties  of  C  are  easy  to  verify. 

(PI):  M  G  C  and  r  real  =>  rM  €  C 
(P2):  M  EC  =>  MT,M*,Re(M)GC 
(P3):  nonsingular  M  G  C  ==>  M-1  G  C 

(P4):  M,N  G  C  =>  M  +  NGC 
(P5):  M,NgC  =>  M0NeC 
(P6):  M,N  G  C  =>  MNMgC. 

Notice  that  MN  need  not  be  a  matrix  in  the  class  C  even  though  M  and  N  are  members 
of  C. 

Recall  that  0  denotes  a  point  between  the  p-th  sensor  and  the  ( p  +  l)-st  sensor  of  the 
array  such  that  (10)  hold  for  some  1  <  p  <  m  -  1  and  0  <  r  <  1.  Clearly,  every  point  on  the 
line  segment  connecting  the  first  and  last  sensors  is  of  the  form  of  O  f°r  some  1  <  p  <  m  —  1 
and  some  0  <  r  <  1. 

We  notice  that 

m  T"X origin  at  O) 

Ao(origin  at  O)  A0(origin  at  0)  =  _  *  (A.2) 

Z)(origin  at  O)  m 

A-l 


where 


m  —  i 

^(origin  at  O)  =  E  el(fc-(p_1+r))Aw 

k=0 

=  .  ei(=±i  — (p+r))A« 


and  hence 


(sin(mAul)\2 

sjn(Zjf  J 

which  is  independent  of  coordinate  reference. 

To  further  simplify  notation  we  write 

£.  (origin  at  O)  =  £(*  ~  (l>  -  ‘  +  O^-O-'+'VW, 

fc=0 

(origin  at  Q)  =  S>  -  (p  -  1  + 


(A.3) 


53fc(origin  at  O)  =  (E(ori6in  at  O))  . 

(or^n  ^  O)  =  (E*(origin  at  O))  • 

We  are  now  ready  to  compute  H(origin  at  O)-  From  (16)  we  obtain 

D2(origin  at  O)  Do(origin  at  O) 

ZT=o(k  -  (P  ~  1  +  r))2  E^(origmatO) 
Zk*  (origin  at  Q)  EETo  (*  ~  (P  ~  1  +  r))2 


where 


m  — 1  i 

(fc  -  (p  -  1  +  r))2  =  -m(m  -  1)(2 m  -  1)  +  m(p  -  1  +  r)(m  -  p  -  r). 

fc=0 

Next,  from  (13)  and  (16)  we  have 

Do(origin  at  O)  Ao(origin  at  0)  =  -iMo(origin  at  0)> 
A2(origin  at  Q)  Do(origin  at  Q)  =  iMo(origin  at  Q), 


(A.4) 


(A.5) 


where 


M0(origin  at  Q)  = 


m(m  +  l-2(p  +  r))/2  £*(ori gin  at  0) 

^(origin  at  Q)  m(m  +  1  -  2(p  +  r))/2 


(A.7) 


EYom  (A.2),  (A.4),  and  (A.7)  we  find  that 


AJKorigin  at  O)  Ao(origin  at  O)  £ 

D5(origin  at  0)  Do(origin  at  O)  €  C, 

Mo(origin  at  O)  €  C. 

By  the  generality  of  the  point  O,  the  matrix  AqAo,  D^Do,  and  Mo  belong  to  the  class  C 
regardless  of  choice  of  coordinate  origin. 

Since 

h  =  dSPaoD°  =  DSDo  -  DSAo(ASAo)-1  ASD0 

=  D5D0  -  Mo(ASA0)-1M0,  (A.8) 

it  follows  immediately  that  H  is  also  an  element  of  C  regardless  of  choice  of  coordinate 
origin. 

Furthermore,  being  Hermitian,  H  is  positive  semidefinite.  That  is, 

Hu  >  0;  H2n  -  |H12|2  >  0.  (A.9) 

We  claim  Hu  ^  0.  To  see  this  we  argue  as  follows.  If  Hu=0,  then  by  the  second  part  of 
(A.9)  we  also  have  Hi2=0  and  thus 

0  =  H  =  D*P^D  =  D*  (Pa)*Pad-  (a-10) 

Let  dj,  and  d2  be  the  two  columns  of  H,  then  (A. 10)  gives 

d;(ri)'Pid.  =  o  =  d;(pi)-pid2, 

•)r(Pi)-Pid2  =  o  =  d;(pi)-pid1, 

which  then  implies  di,d2  G  span(A).  If  r  =  0  (resp.,  r  =  1),  then  the  coordinate  reference  is 
the  p-th  (resp.,  (p+  1) — st)  sensor  element.  This,  however,  is  impossible  since  both  column 
vectors  di  and  d2  have  their  p-th  (resp.,  (p  +  l)-st)  elements  equal  to  zero  while  the  p-th 
(resp.,  (p+  1) — st)  elements  of  a*  and  a2  are  both  equal  to  one.  Thus,  we  can  suppose  that 
0  <  r  <  1.  It  is  clear  that  neither  dj,  j  =  1,2,  is  a  (complex)  multiple  of  ai  or  a2.  Hence,  if 
dj  €  span(A)  then  there  are  two  non-zero  complex  numbers  a  and  /?  such  that 

i(m  -  (p  +  r)  -  j)  =  a  +  /3e^m~^p+r^~^Au/,  0  <  j  <  m  -  1. 


A-3 


This,  however,  is  unattainable  for  m  >  3.  Thus,  Hn/0,  and  hence  Hu  >  0. 

Finally,  using  (A.8)  with  (A.2),  (A.4),  and  (A. 7),  we  obtain  (29)  and  (30)  immediately. 
This  concludes  the  proof  of  Lemma  3. 

For  the  case  of  equipowered  signals  mentioned  at  the  beginning  of  the  appendix,  we  note 
that  SeC  for  tti  =  x2.  From  (2),  (4),  (9)  (for  0  <  |p|  <  1)  and  (6)  (for  |p|  =  1)  we  have 
(A.l)  readily. 
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APPENDIX  B 


A  Nontrivial  Example 

We  shall  consider  the  situation  where  the  electrical  angular  separation  Aw  is  an  integer 
multiple  of  (standard  or  first  null)  beamwidth,  denoted  by  BW,  with  one  BW  equal  to 
2 jt/to.  That  is 

Au  =  l  BW  =  t  . 

Since  0  <  Aw  <  7r  it  follows  that  l  =  1,2, [yj  where  [yj  denotes  the  largest  integer  less 
than  or  equal  to  y. 

The  interesting  result  we  derive  here  can  be  summarized  as  follows. 

Proposition  Append.B  Suppose  m  >  4  and  Au  =  l  BW,  for  t  =  1,2,...,  [yj.  Then, 
for  all  choices  of  coordinate  reference  and  any  signal  correlation  matrix  S,  the  variance 
[CRB]  -  ,  i  =  1,2,  is  finite. 

For  m  =  3  and  Au  —  1  BW  =  the  variance  [CRB]  - ,  t  =  1,2,  is  arbitrarily  large 
(oo)  if  and  only  if  the  two  signals  are  correlated  with  appropriate  phase  difference  0  which 
depends  on  the  choice  of  coordinate  reference. 

Proof  of  this  proposition  is  carried  out  by  the  following  two  lemmas.  It  is  interesting  to 
note  that,  when  m  —  3,  the  electrical  angular  separation  Au  =  1  BW  =  y  is  neither  small 
in  absolute  value  (y)  nor  in  relative  value  (1  BW).  Physically,  for  interelement  spacing 
equal  to  half  of  wavelength,  two  plane  waves  with  angles  of  arrival  (with  respect  to  the 
normal  to  the  array),  say,  <f> i  =  0°  and  such  that  sin(<fo)  =  §  (say,  fa  ~  42°)  give  rise  to 
this  electrical  angular  separation. 

First,  using  explicit  expressions  (29)  and  (30)  for  elements  of  H  given  in  Lemma  3,  we 
have  immediately 

Hi, (origin  at  O)  =  J  ~  ~T^j)  ’  <B  I> 

=  (-1)<mcos(”) 


Hi2(origin  at  Q) 


(B.2) 


hence 


|Hi2(origin  at  0)1  = 


m  cos 


(B.3) 


2  sin2  (£) 

By  the  generality  of  Q,  we  conclude  from  (B.l)  and  (B.3)  that  Hu  and  |Hi2|  are  indepen¬ 
dent  of  the  choice  of  coordinate  reference  when  A w  =  l  BW. 

Since  the  general  coordinate  reference  is  used  in  this  appendix,  to  facilitate  the  discussion 
we  shall  not  restrict  the  temporal  phase  difference  to  any  interval;  however,  in  all  cases,  by 
Lemma  2,  we  know  that  [CRB(0  +  7r)](i  =  [CRB(0)]u  ,  i  =  1,2. 

If  m  =  3  then  £  =  1.  It  follows  immediately 


Hu  =  1  =  |H12|, 

2tt 

arg(Hi2(origin  at  Q))  =  *  +  —  (2  -  (p  +  r)) . 


By  appealing  to  (45),  we  obtain 


Lemma  Append.B.l  Suppose  m  =  3  and  Aw  =  1  BW  (=  y1).  Then,  for  any  choice  of 
coordinate  reference,  the  matrix  H  is  singular.  Consequently,  for  m  —  3  and  Aw  =  1  BW, 
the  variance  [CRB]it ,  i  —  1,2,  is  arbitrary  large  (oo)  if  and  only  if 

\p\  —  1  and  0  =  ~(1  -  2(p  +  r))  -f  kv,  (B.4) 

where  k  is  any  integer,  \p\  and  0  as  in  (26),  and  p  and  r  defined  the  coordinate  reference  0 
as  in  (10).  In  particular,  if  the  center  of  the  array  is  chosen  as  the  coordinate  reference, 
then  9  =  kit  in  (B.4). 


We  next  consider  the  case  where  m  >  4.  For  ease  of  notation  we  shall  define 


H 


ti 


then,  by  (B.l)  and  (B.3), 


f(m;£)  = 


2  cos  ^ 

m) 

1 

(  ms  —  1  ^ 
V  3  J 

|  sin2 

(&)-' 

Suppose  m  >  4  is  fixed.  Then  for  1  <  l  <  it  <  [yj,  we  see  that 


sin(w) >  sin(!) <>o)’ 


(B.5) 


B-2 


and 


(°<)  co.(£)  <  c <*(£). 

Therefon,,  fUr  a  fixed  m  >  4, 

/(m;fc)  <  for  1  <  t  <  k  <  ^yj  .  (B.6) 

We  shall  prove 

f(m-,l)  <1  for  m  >  4  and  1  <  l  <  ^yj  .  (B.7) 

To  this  end,  in  view  of  (B.6),  we  need  only  show 

/(m;  1)  <  1  for  m  >  4.  (B.8) 


From  (B.5),  it  is  clear  that  /(m;  1)  <  1  is  equivalent  to 

1  +  2c°s  (i)  <  sin2  (|)  for  m  >  4. 

We  list  the  coresponding  values  of  the  LHS  and  the  RHS  of  (B.9)  for  4 


m 

l+2cos(%) 

(n^l)sin2(X) 

4 

5 

6 

j  7 

1  +  v/2  »  2.4142 

«  2.6180 

1  +  y/3  w  2.7321 

«  2.8019 

5 

2 

»  2.7639 

f§  w  2.9166 

ss  3.0121 

Observe  that  (B.9)  holds  for  4  <  m  <  7.  Also  note  that  the  LHS  of  (B.9)  is  bounded  by 
3  and  the  RHS  of  (B.9)  is  greater  than  3  for  m  =  7.  Thus,  to  show  (B.9)  it  suffices  to  show 
the  RHS  of  (B.9)  is  an  increasing  function  of  m.  To  do  this,  we  define 

S(x)  =  —3-^  sin2  (y)  ,  *  >  4. 


(B.9) 

<  m  <  7. 


Then 


</'(*)  = 

Using  the  inequality  (reference  [19],  page  75) 


cos(t)  < 


<  8in(0 


0  <  t  <  ff, 
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we  write 


x  . 
—  sin 

T 


>  cos  (! )  -  (^)  “s  (I) 

=  ?cos(?) 

>  0,  for  x  >  4. 


This  proves  g'(x)  >  0  and  thus  g(x )  is  an  increasing  function  of  x,  for  x  >  4.  It  follows  that 
(B.9)  holds  for  all  m  >  4  which  proves  (B.8)  which,  in  turn,  implies  (B.7). 

We  summarize  the  discussion  in 


Lemma  Appd.B.2  Suppose  m  >  4  and  Ao>  =  l  BW,  for  t  =  1,2,...,  [yj.  Then  the 
matrix  H  is  nonsingular  for  all  choices  of  coordinate  reference.  Consequently,  the  variance 

[CRB],-,- ,  i  =  1,2,  is  finite  for  any  choice  of  coordinate  reference  and  any  signal  correlation 
matrix  S. 
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Computing  [CMu(0)]a  ,  i  =  1,2 

In  this  appendix,  we  snail  compute  [CMt/(0)]i;  ,  t  =  1,2,  for  the  center  of  the  array  as 
the  coordinate  reference. 

Recall  that,  for  0  <  |p|  <  1, 

CMu  =  y(H0I)_1  Re(H0KT)  (H0I)_1’ 

where 

K  =  S'1  +  er2S-1  (A*  A)-1  S”1. 

Using 

<H0I>"  -  (jst)  *■ 

together  with  the  fact  that  H  is  real  when  the  coordinate  reference  is  chosen  at  the  center 
of  the  array,  we  obtain 

cmu  =  ^iH0Re(Kr). 

Therefore,  for  i  =  1,2, 

[CmuWI.  =  ^  Kk%  =  ^  l«*(K)h  =  ^Re(K„),  (C.1) 

where,  by  a  straightforward  computation, 


Re(K„) 


1 

(1  ~  I P\2) 


i,j  e  {1,2} 

1  ;(C.2) 

j  /  *• 
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